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ABSTRACT 

Despite extensive study, the mechanisms by which Be star disks acquire high 
CN ' densities and angular momentum while displaying variability on many time scales 

^ ■ are still far from clear. In this paper, we discuss how magnetic torquing may 

CN ■ help explain disk formation with the observed quasi-Keplerian (as opposed to 

expanding) velocity structure and their variability. We focus on the effects of 
the rapid rotation of Be stars, considering the regime where centrifugal forces 
provide the dominant radial support of the disk material. 

Using a kinematic description of the angular velocity, Vfp{r), in the disk and 
^ ■ a parametric model of an aligned field with a strength B{r) we develop analytic 

^ ■ expressions for the disk properties that allow us to estimate the stellar surface 

' field strength necessary to create such a disk for a range of stars on the main- 

sequence. The fields required to form a disk are compared with the bounds 
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CN . previously derived from photospheric limiting conditions. The model explains 

^ . why disks are most common for main-sequence stars at about spectral class B2 

Q_|. V. The earlier type stars with very fast and high density winds would require 

Q . unacceptably strong surface fields (> 10'^ Gauss) to form torqued disks, while 

^ . the late B stars (with their low mass loss rates) tend to form disks that produce 

. only small fluxes in the dominant Be diagnostics. For stars at B2 V the average 

surface field required is about 300 Gauss. The predicted disks provide an intrin- 
^ , sic polarization and a flux at Ha comparable to observations. The radial extent 

of our dense quasi-Keplerian disks is compatible with typical estimates. We also 
discuss whether the effect on field containment of the time dependent accumu- 
lation of matter in the flux tubes/disk can help explain some of the observed 
variability of Be star disks. 

Subject headings: stars: emission-line. Be - stars: magnetic fields - stars: rotation 
- stars: winds, outflows - circumstellar matter - polarization 
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1. Introduction 

The presence of disks around rotating stars has long been reahzed from observations, 
but has never been explained in a satisfactory way. Particularly pronounced are the disks 
around the emission hne Be stars. These can be detected because of the strong double-peaked 
emission lines of Ha, large IR continuum excesses, and intrinsic polarizations. Because 
of these interesting and easily detectable properties. Be stars and their disks have been 
intensively studied for many years (see reviews in Smith, Henrichs, & Fabregat (2000) and 
Yudin (2001)). Nevertheless, the origins of their high densities, angular momentum, spatial 
structure, and variability remain a major puzzle. In addition to the dense equatorial disks, 
it is well known that Be stars have substantial winds, roughly isotropic, with typical M ~ 
XO^(io^^)M0 yr"^ and terminal velocities near v^o — 1000 kms~^ ( Marlborough & Peters 
1986; Grady, Bjorkman, & Snow, 1987; Prinja 1989) 

Hot stars exhibit rotational line broadening indicative of surface equatorial rotation 
speeds which are high, but less than about 80 percent the centrifugal (Keplerian) limit. 
Thus rotation alone is not sufficient to centrifuge material off the equator to form a disk. 
Owing to the high luminosity of the stars, wind material leaves the general stellar surface 
due to radiative forces on the line opacity. Line-driving, however, only works well where 
the wind material is not too optically thick, and has a sufficient velocity gradient for the 
absorbing atomic lines to sweep through the stellar continuum spectrum. Observational 
evidence suggests that Be star equatorial disks are too dense and too slowly expanding for 
line driving alone to support them against gravity. 

One possible mechanism proposed to create disks is the Wind Compressed Disk (WCD) 
model of Bjorkman & Cassinelli (1993). In this, stellar wind matter accelerating out from an 
intermediate latitude on a rotating star has a trajectory orbital plane which crosses the equa- 
torial plane. The continuous outflows of such matter from the upper and lower hemispheres 
'collide' in the equatorial plane and form a compressed disk with density enhancement by 
a factor of ~ 100, which is on the low side for what is needed observationally. Irrespective 
of other theoretical issues raised concerning the WCD model (Owocki, Cranmer, & Gayley 
1996), the fact is that all the observational indications seem to point toward a disk which 
is mainly moving azimuthally, supported against gravity by centrifugal forces, i.e., a near 
Keplerian disk [vr <^ ^ y/GM/r) rather than one which is mainly outflowing {vr ^ Vrf). 

The near Keplerian and high density disk inferred from observations poses questions 
which we address in this paper: How does the material acquire more angular momentum 
than it possessed when it left the stellar surface? How does such an equatorial disk accu- 
mulate material, given that its density and small radial velocity gradient prohibit radiative 
driving from overcoming gravity in the equatorial plane? Proposed answers to these two 
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questions have included schemes in which matter is ejected from the star selectively in the 
pro-grade direction to boost its rotation up to orbital values. In particular, Hanuschik (1999) 
considered how an outflow of particles with an isotropic Maxwellian distribution could create 
an orbital regime comprising the faster pro-grade particles with the rest falling back. Owocki 
& Cranmer (2001) proposed non-radial pulsations with an anisotropic radiation force favor- 
ing pro-grade gas flow. While these ideas are interesting, the direct centrifuging action that a 
strong enough stellar magnetic field could provide seems to be a more direct way of producing 
the disk. 

The possible role of magnetic fields in Be disk production has been mentioned quite 
often (Friend & MacGregor 1984; Poe & Friend 1986; Ignace, CassincUi, & Bjorkman 1996). 
These dealt with apphcations of magnetic rotator wind theory, which explored the various 
degrees of control of the outfiowing gas by the magnetic fields (Belcher & MacGregor 1976; 
Hartmann & MacGregor 1982; Poe, Friend & GassineUi, 1989; Mestel 1990). The magnetic 
rotator models considered the effects of the field on the azimuthal motion of the gas, the 
terminal speed of the wind, and the equatorial mass loss rate. 

Several papers have discussed aspects of the time-dependent accumulation of stellar 
wind matter channeled by strong magnetic fields. This has been mainly in the context of 
early type Chemically Peculiar (CP) stars with very strong magnetic fields (Havnes & Goertz 
1984; Babel & Montmerle 1997a), but with ideas being extended recently toward applications 
to other hot stars and Be star phenomena (Babel & Montmerle 1997b; Donati et al. 2001). 
Specifically, Havnes and Goertz (1984) discussed the effect of the growing 'magneto-spheric' 
density and mechanisms which might disrupt the field, releasing energy and producing X- 
rays. Babel and Montmerle (1997a) applied these ideas to develop a magnetically confined 
wind-shocked model for X-rays from the Ap star IQ Aur and extended this (Babel and 
Montmerle, 1997b) to the young star 9^ Ori G. Donati et al. (2001) applied the model to /3 Gep 
which has a magnetic field much lower than those on CP stars, and which, they concluded, 
is highly oblique to the axis of slow rotation of this star. All of these papers, however, 
neglected the contribution of rotation to the dynamics of the magnetically channeled wind. 
This approximation cannot be valid for stars of near critical rotation since the rotational 
energy density of the channeled wind, and even more so of the shock compressed disk, may far 
exceed the outflow energy density. Nor can it say anything about how disk material acquires 
the quasi- Keplerian angular momentum values observed. In this paper we therefore review 
the problems and address the other limiting regime where the rotational term is dominant 
(centrifugal force balancing gravity) . We do so to evaluate what field is necessary to provide 
sufficient torquing to produce a dense disk with quasi-Keplerian speeds. 

We make use of the fact that near main-sequence B stars have a mass outflow and the 



-4- 



channeling of that flow can lead to shock compressed disk material. The relatively high 
densities in that disk lead to the strong Ha, excess free-free fluxes and electron scattering 
polarization for which Be stars are known. 

We develop equations to find the field needed both to transfer angular momentum to 
material that is driven up from the star by line driven wind forces, and to redirect that flow 
toward the equatorial plane. This regime of a magnetically directed flow is non-spherical 
and, because there is a build up of matter, is also non-steady. The directed flow results in 
the formation of a disk-like structure the azimuthal motion of which is 'quasi-Keplerian'. 
By this we mean that the velocity distribution is not strictly Keplerian but the motion is 
primarily azimuthal and with speeds of order the Keplerian value. 

We have chosen to use an aligned field of dipole-like shape for the basic magnetic topol- 
ogy, B{r,6), with the r-dependence parameterized. As for the gas, it will for the most part 
be assumed to be confined by the field. Our main interest is not in the flow trajectories 
v(r, ^,0), but rather in the azimuthal velocity, v</,((i), in the equatorial plane at radial dis- 
tance, d. The closed field will channel the gas to a specific radial distance, d{l), that depends 
on the latitude I the field tube had at the stellar surface, the larger I corresponding to larger 
d. Field lines from high latitudes, I ^ 7r/4, near the pole, would intercept the equator at very 
large distances, but before reaching those distances the wind energy density will exceed that 
of the field and no longer be magnetically dominated, since even a dipole field declines as d~^. 
Thus, there will be a broad, open, polar wind that fills the envelope volume far beyond the 
disk and above the closed lines of force that arc channeling lower latitude material toward 
the disk. For our study, we consider both the strictly closed field regions and the transition 
to the open ones as illustrated in Figure 1. 

The simplest picture of the situation would be an abrupt "switch" model, in which 
the field is assumed to be fully capable of causing solid body rotation, v^{d) oc d, out 
to an Alfven distance, at which the field energy density has decreased to that in the gas 
flow. Beyond that transition radius the gas would be dominant and in particular its angular 
velocity would satisfy conservation of angular momentum, v^ocl/d. In this representation, 
there is a discontinuous shift in i'^(d) (from oc d to t;^ oc d~^). Although this picture may 
be useful for understanding some basic concepts involved in the torquing of disks, we know 
from magnetic rotator stellar wind theory that no such abrupt shift occurs at the Alfven 
point in either or the angular momentum. 

To allow for a more realistic behavior of the angular speed and angular momentum 

transfer, we introduce a "smooth transition" model, in which a continuous parametric func- 
tion, Vif,{d) is specified for the angular velocity as a function of distance in the equatorial 
plane. This function reaches a maximum, but the angular momentum continues to increase 
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with increasing radius by a factor of about two, transmitted by the bent field hnes. We are 

interested in cases in which Vff,{d) reaches at least to the Keplerian speed vk — ■\/GM/d. 
The field lines will not be considered drawn open until at least a distance, (iouter, where 
the Keplerian ratio, i.e. v^{d) /vK{d) is at its maximal value. The maximal Keplerian ratio 
radius is identified with the Alfven radius and, from that radius outward, we assume that 
the gas motion then tends to the field-free flow limit. Thus, at large o?, it obeys conservation 
of angular momentum, and the particles remain in their orbital planes defined at the point 
where they became free from the magnetic field forces. 

The subsequent flow trajectories in the far field region become similar to those in wind 
compression theory, either in the form of the WCD model or the Wind Compressed Zone 
(WCZ) model of Ignace, CassineUi, & Bjorkman (1996). The WCZ model was developed to 
describe the structure of a wind which has a rotation that is too slow to form a shock com- 
pressed disk, but nonetheless can produce density enhancements in the low latitude regions. 
The enhancements in density (by a factor of 10 or so) have an effect on the observational 
properties of a wide range of stars. The field lines are drawn out in the far field regime as 
described by Ignace, Cassinelli, & Bjorkman (1998). Neither the WCD nor the WCZ models 
can explain the quasi-Keplerian disks of Be stars. 

The relation to WCD theory of the current model for a shock compressed disk is as 
follows. The WCD model is "kinematic" in that the radial velocity structure of the outfiow 
is chosen to obey the well-known beta law distribution, instead of being computed from the 
force of radiation on line opacities. In WCD theory, the velocity component is assumed 
to be determined by conservation of angular momentum, and thus f ^ decreases rapidly with 
radius and does not attain the quasi-Keplerian angular speeds inferred from observations. 
In the current paper, we again use a kinematic approach, but now specify an azimuthal 
velocity structure, v^{d), that assures quasi-Keplerian velocities are reached for a sufficiently 
strong magnetic field. This kinematic approach allows us to derive constraints on the surface 
field. Bo, and the surface spin rate parameter, <So, that are needed for a disk to form. Mass 
and magnetic flux conservation requirements provide information regarding the disk density 
distribution, the radial extent, and the timescale for variability. Where needed we also use 
the wind beta velocity law to provide an estimate of the flow velocities driven by the star's 
radiation field. The derived structure for stars with "torqued disks" allows us to estimate 
observational properties such as Ha line strengths, intrinsic polarizations, and continuum 
fiuxes at a variety of wavelengths useful for comparisons of the model with observations. 

A major goal is to derive lower limits on the field strengths required for the existence 
of centrifugally-supported disks produced by an afigned dipole-like field on a Be star. As a 
first rough estimate, showing that the fields required to channel a fiow are not unreasonable. 
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we consider magnetic rotator wind theory. In that theory, as explained in Ch 9 in Lamers 
& CassineUi (1999) (hereafter L&C), there is a "primary mass loss mechanism" (in our case 
the line driven wind theory), which sets the stellar mass loss rate, and the terminal velocity. 
An amplification of both M and v^o can be produced by the co-rotating magnetic field. 
A transition occurs from the case in which the field acceleration effect is unimportant (a 
"slow magnetic rotator") to the case in which the field produces a strong effect (a "fast 
magnetic rotator", or FMR). This transition occurs when the Michel velocity, Vm, (which is 
very close to the terminal velocity in FMR theory) is as large as the terminal velocity owing 
to the primary mass loss mechanism. The Michel velocity is related to the Poynting flux 
which makes the enhanced outflow energetically possible. As a starting point, we consider 
it reasonable to assume that the field necessary to dominate the fiow and channel it to the 
equatorial disk should be at least comparable to the transition field strength from slow to 
fast magnetic rotator theory for a particular star. The Michel Velocity, Vm-, is given by 
quantities defined at the surface of the star; the surface field, (-Bo), rotation rate (r2o), and 
mass loss rate, M. Solving for Bo in the equation for Vm gives a minimal field based on 
FMR considerations. 
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where it! is the stellar radius. We set RQ,o equal to the surface equatorial speed Vo, which 
is expressed as the equatorial Keplerian speed, i'ft'(R) times a fraction So. We also set 

Vm — ~ (where ve is the escape speed at the stellar surface). For numerical values, we 
express the velocity by the ratio vg = v^o/ (10^ cm s^^); the stellar mass and radius ratios by 
Mio = M/(10Mq), Ri2 = i?/(10^2 cni); and the mass loss rate by M_9 = M/{IQ-^Mq yr"^). 

A roughly similar value for the magnetic field, as in Equation (3), follows from requiring 
that the energy density of the surface field Sq^/Stt should exceed that of the gas pov1/2. The 
derived field of order 10 Gauss (for So— 1) is a very modest one. However, as we shall see, it is 
a lower bound on Bo because it effectively uses the density in the wind, which is much below 

(i.e. ~ 10^'^-lO^^x) the densities that arc produced in the centrifugally supported disk, so 

1 

that the field (oc p2) torquing the disk is 30 to 100 times larger than that in Equation (3). 
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2. Schematic Description of the Model 

Let us start by considering the case of a non-rotating star with a strong surface magnetic 
field, which we take to be of dipole-hke form (but of more a general r-dependence) , and with 
a uniform wind mass flux, Tu,o = M/(47ri?^), from its surface. If we consider starting from 
a vacuum environment and initiate an outflow, gas will be channeled along the field lines 
(assuming |p v^j <S B^/Stt), and along with matter from the opposite hemisphere, tend to 
fill the fiux tubes to ever greater distances with wind compressed gas, creating an extended 
dense disk. This is the non-rotating case considered by Babel & Montmerle (1997a). The 
flux tubes would fill toward pressure balance. However, the resulting decrease in the velocity 
gradient would probably quench radiative driving, tending to cause the matter to fall back, 
although Babel & Montmerle (1997a) suggest otherwise. We suspect that the Rayleigh 
Taylor instability might be important and prevent a wind flow from supporting a relatively 
dense envelope against gravity. 

When we introduce rotation, the magnetic field exerts an azimuthal force on the gas as 
well as acting to constrain its radial fiow path to be along the fiow tubes. As is illustrated in 
Figure 1, the latitude / at which the fiow originates will determine the azimuthal velocity of 
matter as it is channeled into the equatorial disk. Matter originating in some latitude range 
li < I < I2 will tend to reach the equator at dinner < d < douter with moderate speeds v\\ 
parallel to the field lines in the corotating frame. As a result of the infiux of matter, there 
will be a build up of mass in these sectors of the disk. Although compressed to high densities 
by the shock fronts of the incident fiows, the material in the disk will be supported against 
infall mostly by centrifugal forces. 

We discuss below the parameters of the various resulting regions in the equatorial plane, 
and also the timescale on which they change. We find the latter to be long compared to 
typical rotation periods. So, for the purposes of this paper, we assume we have reached 
a nearly steady stage of moderately dense disk creation and consider the effect of stellar 
rotation. 

Most of our interest is in the collision of winds in the equatorial plane, which involves 
material that is near and above the local Keplerian speed. Therefore, we find it convenient 
to express equations in terms of the dimensionless equatorial distance variable x = d/R, 
and the variable S{x) = v^{x)/vk{x)^ which is the ratio of the azimuthal speed to the local 
Keplerian speed vk{x). A fundamental parameter for a specific star is the "stellar spin rate 
parameter" <So, which is the ratio Vo/vk{R) of the stellar rotation speed to the Keplerian 
velocity at the stellar equator; thus So— <S(at x — 1). In principle. So can have values in the 
range < <So < 1. However, for most Be stars the range is 0.6 < <So < 0.8. We assume for 
simplicity that the star is spherical and not distorted by the rotation. 
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3. Existence and Extent of Magnetically Centrifuged Disk 

We begin by finding a form for v^{d) tliat is consistent witli known magnetically cen- 
trifuged wind systems and provides a kinematic distribution that we can use at all distances. 
Although the topology of the field in the rotating dipole case is very different, there are 
properties of the azimuthal velocity that should carry over from magnetic rotator theory 
(L&C Chapter 9). At small radii, the azimuthal speed behaves as in the solid body case, 
i.e. V(j) oc d, while farther out it satisfies the conservation of angular momentum relation 
Vrf, OC 1/d. Now let us find a kinematic expression for V(j„ that is valid both near and far from 
the star and also provides a continuous form that asymptotically matches both behaviors. 



3.1. The Description of the Angular Velocity, v^{x), in the Disk 

At small distances from the star where the field is strong, we expect the azimuthal speed 
of disk matter to be controlled by centrifuging, i.e., to follow the co-rotation distribution 

V^CriVo X {X> 1) (4) 

Farther out, the angular velocity in the region well beyond the Alfven distance should follow 
the conservation of angular momentum expression = (v0,peak dpeak)/d, where v^p^peak is the 
maximum V(j){d), which occurs at d — dpeak- We express dpeak in terms of the dimensionless 
parameter, a, where 

a = dpeak/ R = a^peak (5) 

and thus f</,,peak ~ VoCi. At large x we will have 

2 

V a 

v^{x) ~ — (x > 1, a; > a^peak) (6) 

At intermediate equatorial distances, the exact form of Vfj,{x) depends on details of how 
the field bends in reaction to the inertia of the matter. However, a reasonable form for v^{x) 
parameterized by a, and which has the limiting behaviors (4) and (6), is 

where the numerator normalizes v^{x) so that v^{x — 1) = Vo- This analytic form proves 
convenient for exploring the properties of the disk driving regime. Although Equation (7) is 
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purely empirical, it is designed to reflect the actual behavior of matter torqued in a magnetic 
field. 

Figure 2 shows v^{x) from Equation (7) for various a values, as functions of x. Speeds 
are shown relative to Vo, so the relative height of v^{x) and vxix) depend on So- 

While it is interesting to see the dependence of the adopted v^{x) distribution on a, 
we will see later that a is not a fundamental stellar parameter comparable to the rotation 
quantity <So. In fact, for a star with a given M and Voo, the location x — a oi the peak of 
the curve itself depends on So as well as on the field Bo at the surface of the star. Only after 
we explore the consequences of assuming Equation (7) for v</,(x), and particular forms for 
the disk field and density distributions, will we be able to identify the various independent 
stellar parameters for the Be star disk problem. 

Two other curves of interest are also shown in Figure 2: the Keplerian circular velocity, 
vk{x) and the escape speed ve{x), where 



1 , , IGM Vo , . 

M^) = ^v4x) = ^^^^ (8) 

We need to compare the angular speed, v^f, with the local Keplerian speed, so it is 
convenient to use the Keplerian ratio {S{x) — v^{x)/vk{x)) introduced previously. Using 
Equation (7) we get 



S{x) ^ v^{x)x^''' ^ {a^ + 
So Vo + 

This reaches a maximum at a; = \/3a. Figure 3 shows S{x) versus x, for various a values. 
Of particular interest is the value of a for which S just tangcntially reaches unity, i.e. for 
which there is just a single point at which the angular speed reaches the Keplerian speed. 
This "Threshold value" of a is called "ajh", and it occurs at xjh = VSajh. For a given star, 
ajh will determine the minimal magnetic field that is required to produce a torqued disk, 
i.e. one with material having angular velocities at least as large as Keplerian speed. 

Figures 2 and 3 reveal a number of regimes of behavior according to the value of a, 

1. For the 'weak' field regime, curves of type /, with («< ajh) never produce a Vip{x) value 
as high as vk at any x, so <S < 1 at all a;. The centrifugal force supplied by magnetic 
torquing to material originating at any latitude is not able to overcome gravity. We 
assume that the matter will tend to fall back to the star unless supported by other 
forces such as radiation pressure or ram pressure of the inflowing accumulating gas. 
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Field lines will already be bent by the inertia of the gas at radial distances close to the 
stellar surface. 

2. For a "moderate" field regime, for which ajh ^ ct ^ 3 (depending on So) is of most 
interest in producing quasi-Keplerian disks. The surface fields at intermediate I are 
strong enough to centrifuge gas to equal (critical curve A) or exceed (curves of type //) 
the Keplerian circular speed vk{x). This range in a will result in a quasi-Keplerian 
disk over a range dinner to ajoutcr- The distance Xinncr is the first point at which the 
material has a speed greater than the Keplerian speed, and we choose Xoutcr to be the 
distance at which the Keplerian ratio S is maximal and at which we can consider the 
field to lose its dominance over the gas flow. 

3. For the "large" magnetic field regime (curves of type /// for which a » ajh), the 
magnetic torquing gives the matter sufficient v^f, alone to equal escape speed. However, 
the field could still dominate and prevent outflow. Beyond some large distance, outflow 
is assumed to occur with little further increase in angular momentum and the angular 
velocity will vary as x~^. Arbitrarily large a could centrifuge and contain a disk of 
arbitrary radial extent. However, a depends on the surface field for which there are 
bounds as described by Maheswaran & Cassinelli (1988, 1992), so there are in fact 
upper bounds to a as well. 

We should also emphasize that even though the curves of in Figure 2 have the 
appearance of trajectories, they are not matter trajectories, but rather give the azimuthal 
velocities of matter which has arrived at a particular distance x in the equatorial plane, 
having originated from latitude I. 



3.2. Boundaries and Regions of the Disk 

Here we derive expressions for the limiting radii of particular interest to the torqued 
disk problem. The innermost distance of interest, Xinner, is that at which the torqued matter 
first attains the local Keplerian speed. It occurs where S — 1. So by Equation (9), dinner is 
the solution of 

fa^ + iW- 3/2 1 
a^ + x- ~ " S~ ^ ^ 

This non-linear equation can be solved by a Newton iteration, which converges quickly 
when we start with the large a solution to Equation (10), i.e. Xmner — So~'^^^. Figure 4 shows 
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2^inner(Q;, »So) vcrsus CK for various So, (0.5, 0.7, 0.9). 

The next distance of interest is where Vfj, peaks, namely lipeak = a^peak R, which we have 
already seen is given by 



a^peak = tt (11) 

and at which the rotation speed ratio is 

V<p{Xpeak) _ + 1 



Vo 2a 



(12) 



To define an outer radius of the torqued disk we need a criterion to decide where the 
quasi- Keplerian disk torquing effectively ceases. For this we adopt the distance at which 
the ratio S{x) maximizes, i.e. where dS/dx — 0, for a fixed a. Prom differentiation of 
Equation (9) this is found to be 



a^outer = "\/3q; (13) 

The value of at Xouter is given by 



Vo 

which for large a varies as ~ 0.43q;. Thus, both the location and velocity at the peak increase 
roughly linearly with the parameter a. Note that although Xouter is not explicitly dependent 
on So, it depends on So through a{So). 

A key issue is the minimum value of a for which a Keplerian regime exists at all. There 
is a single point at which the angular speed just reaches the Keplerian speed. Setting S — 1 
X — -\/3q; in Equation (9) gives 

"th"^ + 57^ = 53^ (18) 

This minimal or "threshold" value of a, produces a Keplerian value of V(f, at the singular 
distance xjh = V^ctjh- In Figure 5 wc show the numerical solution for aj\^{So) obtained from 
Equation (15). This equation shows that ajh will be less than 1 for values of So larger than 
2/(3^/^) ^ 0.9 (also shown in Figure 5). However, it should be noted that at very large spin 
rates the star will be substantially non-spherical and our whole treatment becomes rather 
approximate. Our results should be regarded as most reliable for values of So that are less 
than 0.9. 
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We observe from Figure 5 that for the So range most relevant to Be stars, an excellent 
fit to the solution is 

OiTh ~ So'^; where K ~ 0.88 (16) 
shown as the dashed line in Figure 5. 

To relate a to absolute scales of B.p, and So, wc will adopt x = Xoutcr = a/S^ as the 
point where B'^/Stt ~ i.e. wc take Xoutcr to be the Alfven distance. Therefore, our 

definition of the Alfven distance leads to the relation 

XAii = a^outer = a/SQ! for a ^ ajh (17) 
Thus, the Alfven distance is the outer boundary of our torqued disk. 

We recognize that there is a degree of arbitrariness with this definition of the Alfven 
distance. However, just as we were guided by wind theory in our choice of the kinematic 
expression for the angular velocity, we can also get some insight from magnetic rotator winds 
regarding an appropriate definition of the Alfven radius. In the case of winds the Alfven 
radius does not correspond to the peak of the f(^(x) distribution, but as seen in L&C, Fig 
9.3, occurs farther out because the field continues to contribute to the speed of the gas. 
The peak of our S curve (at x = \/3a) clearly lies farther out than the maximum in the 
curve (at x — a). It corresponds to an extremum in the torquing of the disks, in that 
no further increase in the Keplerian ratio occurs. However, perhaps the best reasons for 
choosing the maximum in S to be the outer boundary of a disk are seen by considering 
the threshold case in which the S curve just tangentially reaches unity, which occurs for 
a = ajh and at the location, .xjh- By definition of "a threshold distance" we would want 
the inner and outer boundaries to coincide, and the condition dinner = ^^outcr occurs at xjh- 
The tangentially Keplerian point is also an important location because the components of 
the total velocity {vr, vo,v^) are fully definable for the first time at xj^: The radial velocity, 
Vr equals zero at the equatorial extremity of the dipolar field. The polar velocity vq is there 
nulled by the coUision between the oppositely directed streams. The tangential velocity there 
is v^i, — vxixTh), so Vtctaii^Th) — ^'(/.(^^Th)- Furthermore, using shock conditions we will find 
the distance xjh to be the first radial distance in the equator for which our marginal model 
leads to a quasi-Keplerian disk density (po)- Therefore, all the quantities needed for giving 
an equality between the total speed and an Alfven speed, B / ^Anp, are definable at xjh- 
Thus we will be able to use the conditions at Xjh to define a threshold surface field i^o,Th 
required to produce a torqued disk. 

Choosing the peak of S to define the outer boundary of the torqued disk and also the 
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Alfven radius will allow us to predict the dependence of observational fluxes on field strengths 
and intrinsic stellar properties. We will then be able to test the resulting models against 
general Be and hot-star observational requirements. 

We assume that, beyond Xouter(= the gas flow will tend to draw the field outward 

and change from being predominantly an azimuthal flow to dominantly an outflow in the 
form of ficld-frcc WCD or WCZ trajectories. The flow will make the transition to a nearly 
radial wind as the latitude of the point of origin increases. For a > ajh, the 'Keplerian disk' 
regime extends from d = Xinncr R (Equation 10) to d = XoutcrR (Equation 13). Beyond the 
distance .Xoutcr = V^a, the flow is no longer controlled by the flcld and there is a transition 
from a quasi Keplerian to an outflow dominated velocity structure. 

To get numerical values for these various boundaries in x, we need to specify how the 
magnitude of B{x), and p{x) vary with distance, for a given value of So- Since a depends on 
(So, as well as on B and p, we specify all of these quantities in Section 4. To get a rough idea 
at present of the expected disk extent, we take <So=0.7 and assume that B and p are such 
that a — 2. Applying Equations (10) and (13), we find the range of the quasi-Keplerian disk 
to be 1A5R < d < 3.46it!; a plausible range for Be star disks. 



We choose to represent the field in the system as having a topology similar to the dipole 
limit, but to allow for greater generality we use a parameterized rate of decline with distance 
d. Specifically we set 



where Bo = B{R) is the value at the stellar surface. In fact, the disk does not extend down 
to a; = 1 but only to a; = dinner- Beyond x = Xoutcr the fleld lines tend to stretch and break 
open and we do not attempt to describe this here as it is not inside the torqued disk regime. 
As we will show, the field within the disk determines its radial extent. We also need to know 
the shape of the field fines in the region out of the equatorial plane. This is because it is 
the frozen-in tube geometry between the footprint of the field at stellar latitude I, and the 
point x{l) where the field crosses the disk which determines the flow geometry (i.e. area of 
the tubes) through which the wind mass fiux is channeled. 



4. The Distribution of B{x) and p{x) in the Equatorial Disk 



4.1. Distribution of B{d) 





(18) 
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To get the tube geometry, we make the simphfying approximation that the field strength 
(and also the wind mass flux) is uniform over the surface of the star, i.e. is assumed to 
be the fleld at all /, and the fleld lines are normal to the stellar surface. These assumptions 
— as well as constant radial B and constant mass flux - will not be strictly true in a rapidly 
rotating star, but they should be adequate to describe the limited range of latitudes / that 
adds matter to the range dinner to Xouter in the disk. 

With these assumptions regarding the field geometry, the outgoing cross-section of a fiux 
tube of area Ao at the stellar surface has a value ^4 at o? in the equatorial plane determined 
by field line conservation to be 

A{d) = Ao ^ = Ay (19) 

where we have used Equation (18) for the ratio of the fields at the base and at the point of 
intersection with the disk. 



4.2. The Density Structure, poid), of the Disk 

We assume that the density is determined by pressure balance in the post shock region 
where the flows from the upper and lower hemispheres collide. Thus, the density near the 
equatorial plane is determined in much the same way as in Bjorkman & Cassinelli (1993), 
but the material is not flowing outward, in contrast with the WCD model. Also, since the 
material has an angular velocity at or above the Keplerian speed, it is supported against 
gravity by centrifugal forces. Thus the dynamic pressure does not also need to support the 
gas against gravity, which is the case for the models of the slowly rotating chemically peculiar 
stars (Babel & Montmerle 1997a). 

The opposing wind streams provide a ram pressure that must be balanced by the static 
gas pressure in the shock compressed region. The density in the shock compressed disk can 
be significantly larger than that in the incident winds, especially in the near-equator zone 
where the shocked gas has cooled radiatively. We can assume the equatorial disk region to be 
cool 10^ K) because such a temperature can be maintained from a radiative equilibrium 
balance with the incident stellar radiation field. That is, the density is sufficiently high that 
radiative cooling is effective and sufficient to balance the heating by the stellar fight. For a 
Be star, the appropriate temperature range is from 1 to 2 x 10^ K, and with the corresponding 
sound speed cs — {kT/ fXTJiHY^^, where fJ^mn is the mean particle mass. We find by equating 
the disk gas pressure to the wind ram pressure that. 
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Pd{x)(^s = Pw{x)vy,'^{x) = Tm{x)v^{x) (20) 

where pu, and Vw are the density and velocity of the flow at the equatorial end of the magnetic 
tube, and J-'m{x) is the mass flux there. We use the subscript, w (for wind), on both the 
incident velocity and density. We assume that the radiation forces drive a field-channeled 
fiow with a parallel velocity, which is comparable to the wind speed that could have been 
reached at a radial distance d in the line driven wind of a non-magnetic B-star. We take this 
to be given by the beta law expression 



vM) = ^oo (^1 - (21) 

The mass loss rate of the B-star sets the surface mass fiux scale J-m, and as stated above 
it is sufficient for our purposes to assume the mass ffux rate J-'m,o = ^m{x = 1) is the same 
at all latitudes on the star. The mass ffux at x is related to the mass ffux entering at the 
star via J-'m{x) x A{x) = J^u,o x Aq, where A is the area perpendicular to the ffux tubes, as 
given in (19). Thus we obtain, 



^m(-) = ^--^ (22) 



= 5.02 X 10-'^(gm cm-2 s"!)— ^ x'^ (23) 

Using Equations (21) and (22), the disk density at the equator is given by the pressure 
balance equation (20), thus 



2 



M V. 



oo 



= Pc(^)'(l-l/x)^X-^ 



(25) 
(26) 



3.04 X 10-13 g^^^-3^^^(i _ ^i^Y x-^ (27) 



where we let 



Pc = M/{AnR\^) = 5.02 x lO'^^gm cm'^^-^ (28) 



-16- 



be a characteristic stellar mass loss density, and 

Pd,c = Pc{vjcsf = 6.06 X lOVc ) (29) 

be the characteristic density in the disk. The corresponding number densities are given by 

Tie = Pc/(/W"^//) and nD,c = PD,c/(At"^i/)- 

From differentiation of Equation (24), we find that the disk density peaks at x = 1 + 
P/b < 1 + P/3 which lies close to dinner, i-e. 



The disk density in Equation (27) is nearly 4 orders of magnitude larger than that in 
the "wind" incident on the disk 



= 5.02xl0-^^gmcm-3i^^(l-l/a;)-'^x-'' (32) 



Thus the two adjacent densities (post- and pre-shock) given by Equations (27) and (31) 
respectively, differ by about a characteristic Mach number squared (t'co/cg)^, as is typically 
the case for isothermal shocks. In the case of early-type stars such as the Be stars, the Mach 
number factor is quite large. Finally, the large density contrast between the region in and 
out of the shocked disk justifies the neglect of gas pressure compared with the centrifugal 
support, as assumed in our model. 

5. The Magnetic Field for Disk Production in Terms of Stellar Parameters 

5.1. The Alfven Distance and the Minimal Field to Produce a Disk 

We now have explicit parameterizations of the disk field, B{x) (Equation 18), density 
Pt){x) (Equation 24), and azimuthal velocity (Equation 7) as functions of x. Given these, 
we can explicitly find the minimum or threshold surface field, So,Th, needed to torque a 
quasi- Keplerian disk, and the extent of the disk, as functions of the stellar spin rate. So, 
and wind properties, M, /3, and v^o- The existence of unique solutions for Sq^th and for disk 
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inner and outer radii depends on the facts that v^{x) increases with x in the inner region 
and so can exceed vk if B is strong enough, and that i?^(x)/47rpD dechnes monotonically 
with X, so that the torquing effect of the field turns off at some sufficiently large x. It is 
easy to confirm that Equations (18), (24), and (7) satisfy this monotonic decline in energy 
density ratio. 

To express a explicitly in terms of the physical parameters of the star, we use the fact 
that at the rotational Alfven point, a;Aif (=a;outer), 

^ ^^^^^"^ = ^PT)ixAli)Vcl>^{xAlf) (33) 

which, on using Equations (18, 24, and 7) implies that 

Using xah = VSa (Equation (17)), and substituting this in Equation (34) gives us 
where we use the definition 



R 4 1 

7 (^s 



(36) 
(37) 



Equation (35) is our basic equation for relating a, the fundamental quantity in our kinematic 
model, to the stellar field and basic stellar properties contained in F. In Figure 6, we show 
F(q;) from Equation (35) for several h values (3, 6, 9) with /3 = 1, while illustrating a range 
of /3 values (0.5, 1.0, 2.0) for the pure dipole-field case (6 = 3). Note that 



7 



Bo 



Stt 



GMPc 
2R 



(38) 



Therefore, 7 is a measure of the magnetic energy density to gravitational energy density of 
the wind material near the stellar surface, and depends only on So, M/R x Mjv^. Thus 7 
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is essentially the ratio of the surface Alfven speed at the density of the wind to the surface 
Keplcrian speed. On the other hand, So is simply the ratio of surface rotation speed to 
surface Kepler speed, so the right-hand side of Equation (35) (i.e. F) involves only velocity 
ratios at the star while the left-hand side contains only the model parameter a plus the 
wind f3 value, which is known {i.e. ~ 1.0), and b which parameterizes the rate of decline of 
the field and we know b > 3. Note in particular that we can identify, 7, foo/cs and So as 
independent and fundamental parameters, whereas the kinematic description parameter a 
is not fundamental, as o; = a{'j,v^/cs,So)- 

5.2. The Fundamental Stellcir Pctrameters for Disk Production 

We can now translate the analysis of the disk, in §3, in terms of kinematic parameter a 
into the real physical parameters of different stars. 

First of all, we can find the surface field Boju needed to create a quasi-Keplerian disk 
for a star of given So, M, R, M, Voo, and f3. To do so we use the results of Figure 5 (from 
Equations 15 and 16) to obtain ajh{So) for the given So value. Wc can then use ajh in 
Equation (35) (with the chosen values of the parameters b, and P) to obtain FTh('5o) which 
is a convenient intermediate variable encompassing all of the velocity ratios of the system. 

Given Fth from aj^ for the chosen /3 and b and using Equation (28) for pc, we can now 



find 



1 




(39) 



So TTh{So) X 



(40) 



where Y* is determined by stellar properties. 



Using the reference values for M, R, M, and v^o introduced for Equation (3) while 
introducing cs — 1.3 x 10^74^^^, Equation (39) can be expressed numerically as 



1 




(41) 



In principle we should solve Equation (15) to get aTh{So) and hence rTh(»5o), but we 
noted that Equation (16) {a — So~^) yields an excellent empirical fit to ajhi^o)- Using the 
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empirical fit in Equation (36) gives 

06/4+1/2 / q K\f^/'^ 

rTh(5o) = ^-^s.-^^'/'-'\s.-'^ + (42) 

which we show in Figure 7 along with the exact solution of Equation (15). 

Inserting K = 0.88 exphcitly in Equations (39) and (42) for the range of So of greatest 
interest for Be stars we get, 



BoMSo) = 17.5 Gauss x 3V4+i/25^i.88-o.44;,(^^-i.76 + 1) _ 



\/3 / I -Rl2^4 



(43) 



Figure 8 shows -Bo,Th versus So for a B2 V star for various values of b and /?. For the 
dipole case (6 = 3) the threshold surface field for a star rotating at «So=0.7 is about 250 
Gauss. If the field drops off more rapidly with distance from the star, a significantly larger 
surface field is required to produce the disk. Primarily this is because the inner boundary, 
which is identical to xau for Boju, is not at the stellar surface but at 2.43 R, so for the 
larger values of b, a larger surface field is needed to give the same field at xau- The effect of 
changing the wind velocity law index is shown by the curves labeled with various values of 
f3 in Figure 8, and these indicate that the steepness of the wind velocity law is also of some 
importance in determining the field strengths required to produce a disk. For example, in 
the case with So=0.7, a slowly-accelerating flow with /5 = 2 can be channeled to the equator 
by a field only 70% as strong as that needed for the more-quickly-accelerating P — 0.5 case. 



6. Properties of a Torqued Quasi-Keplerian Disk 

6.1. Hot Main-Sequence Stars and the Possibilities for Disk Production 

To see how Bojh depends on stellar parameters we note that M, R, M, Voo, along the 
main sequence are all essentially determined by the spectral type Teff as given in a table from 
Bjorkman & CassineUi (1993). Table 1 gives the factor in Equation (35) which isolates the 
field dependence on stellar properties. Adopting b = 3 and f3 ~ 1 we can then find Bojh{So) 
as a function of Tefj from Equation (35). Table 1 also includes representative values of B^jh 
for 5o=0.5, 0.7, 0.9. 

Figure 9 shows the comparison the minimal field from Equation (43) along with the 



-20- 



bounds on stellar magnetic fields that were derived by Maheswaran & CassineUi (1988, 
1992). The latter bounds were derived from considerations of envelope circulation currents 
and maximal contributions of field pressure to hydrostatic support. We can see in Figure 9 
that very large magnetic fields are required to produce disks around O stars. Such large 
fields are needed because the large wind momenta (Mvoo) of 0-star require large magnetic 
forces to deflect their flows. For the case of the 03 V star the required fleld lies above the 
maximal allowed field for many different stellar spin rates. For the early B stars with a spin 
rate of So — 0.7 the field required is in the range 100 to 1000 Gauss. Such fields are small 
enough to have avoided detection in the past, but are near the field strength recently found 
by Donati et al. (2001) in /3 Cep. For the late B stars, we see that because of their small 
wind momentum rates, rather weak dipole fields of order 30 Gauss could produce a torqued 
disk. In such cases we need to know if the predicted disks would be detectable. 

The properties of the stars in Table 1 lead to specific predictions regarding the disk 
extent and observational diagnostics should vary with stellar parameters. In particular, 
for given values of b, (3, M, R, we expect the radial extent of the quasi-Keplcrian disk to 
increase with increased Bo while decreasing with increased Mvoo, predictions that might be 
observationally testable. We now investigate the radial extent of disks as function of stellar 
magnetic field, and find the number of particles in the disk to find if the model can explain 
the polarimetric properties and disk emission measures, which determine the Ha fluxes, of 
Be stars. 



6.2. Radial Extent 

In Section 3 we obtained expressions in terms of a, So, for the radial extent of the disk, 
2^inner a;outer- This cxtcut is of coursc zcro for Bo < Bojh- For Bo > Bojh, we can obtain 
explicit values for Xinner and Xouter, in terms of Bo, So and other stellar parameters (M, R, 
M, Voo, T), by translation from a to F using (38). Speciflcally, for any given set of stellar 
parameters Y^, we can calculate F and then obtain the corresponding a (for given (3, b) from 
Figure 6, or by numerical inversion of (38). Then we have Xouter = xwf — V^a. Here we 
illustrate results by adopting b — 3, (3 — 1, and deriving dinner and Xouter as functions of So 
and the ratio Bo/Y^ 

The results are shown in Figure 10. Figure 10a shows that results for three values of 

So- On this panel, it is easy to see that So sets the threshold for the field needed to produce 
a disk. This figure also shows that for values of B /Y^ slightly above the threshold field, 
the disk extent increases sharply, and in a seemingly homologous way vs. So- To check the 
apparent homology we show in Figure 10b the ratio of dinner and ajouter expressed in terms of 
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threshold radii xju versus the ratio of B/Bojh- 

As for the inner boundary, the position Xi^ner denotes the position in the disk where the 
material first reaches the Keplerian speed (i.e. S — 1). For the minimum surface field which 
will accelerate material to Keplerian speed (Bojh), the disk material reaches Keplerian speed 
at only one point, a^inner = stouter = VSa. The other point on the kinematic Vcj, curve in Figure 
2 where the speed is Keplerian is of no interest because it is beyond the torqucd disk outer 
boundary, where the field is no longer dominant in determining the disk's ^^(x) structure. 

Note in Figure 10 that, for values of Bo large enough to enforce nearly sohd body 
rotation, the location of the inner boundary is nearly independent of the magnetic field, as it 
is primarily set by the value of So as mentioned earlier. In regards to the outer boundary, an 
expression for Bo/Bojh can be derived from the ratio of two values of in Equation (36), 
evaluated at a = xau/V^ and at a = xj^/ \/3, and we find that the outer radius varies with 
Bo/Bojh, roughly as ^Aif/a^jh ^ (-Bo/-Bo,Th)^^^- The net effect is that the outer boundary 
and the extent of the disk increases as Bo increases, but the location of the inner boundary 
does not change much. 

From an observational point of view these considerations only concern the magnetically 
torqued regions while the flow in the far regions could also show enhanced densities because 
of the WCD effects as schematically illustrated in Figure 1. 

The typical range of the disk model boundaries here is of the same order with Xoutcr ~ 3 
to 10, as discussed in the empirical models of Be stars assuming Keplerian disks and based on 
IRAS observations (Waters 1986), and Ha line observations by Cote, Waters, & Marlborough 
(1996). As the latter authors note, the observations do not rule out disks of greater radial 
extent, since the tenuous outer parts contribute little to the polarization or emission measure. 
Indeed some observations such as the radio observations of Taylor et al. (1990) show that 
high density material can reach distances of order 100 R. The bulk of this extent might be 
a tenuous outflow WCD/WCZ region as opposed to a quasi- Keplerian one. 

Now we have the physical boundaries of a Magnetically Torqued Disk and expressions 
for the density distribution in the disk so it is possible to test if the predicted disk satisfies 
basic observational limits, especially in regards to the classic Be star diagnostics provided 
by polarization and Ho; observations. 
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7. Compcirisons with Observations 

7.1. Continuum Scattering Polarization 

One important diagnostic of the Be disk properties is the polarization they produce by 
electron scattering. Using the Brown & McLean (1977) general theory of (single) scattering 
polarization, one finds that for a disk of approximately uniform thickness H and a density 
distribution given by Equation (34) for Xmner < 2; < Xouter, the resulting polarization is 



3 aT /"^outer 

p = — — / ne{x)V{x)Hdx/x (44) 

32 TT /t- 

inner 

Here (Jt is the Thomson cross section, rig is the electron number density (= p-o/f^i^H), 
and we ignore absorption and assume the disk is fully ionized. Here we have included the 
finite source depolarization factor, V of Cassinelli, Nordsieck, & Murison (1987), and Brown, 
Carlaw, & Cassinelli (1991), and assumed an edge-on aspect angle. Equation (44) yields, for 
ne{x) given by Equation (24) 



P = ^^^D,c/p (45) 



where the integral 

/l/-^inner 
y^-\l-yf+y^dy (46) 
-/•^ outer 

can be expressed in terms of incomplete Beta functions. 

For Xouter > 1, /9 = 0.5, and 6 = 3, the integral Ip fa l/(2xinner^) - l/(3xouter^) ~ 0.1 
for typical dinner- Using Equation (24) for the disk density in Equation (45), the resulting 
numerical value for the polarization is given by: 



'nB,c\ „ H 
— 77 R12 X — 
.IQii/ R 

Using the theoretical disk number density from Equation (24) 



p(%)~ 0.5(^)i?i2X- (47) 



n{x) = rzD.c ^-''(1 - 1/xf = 1.8 x lO^W^^ ^"^^^ ^^^^^ ^^""^ ^ (48) 

then Equation (47) produces p ~ 1% for the scale values M_9, i?i2, fs, T4 all = 1. Polariza- 
tions of about 1% are typical for the intrinsic polarizations of Be stars (McLean & Brown 
1978). 
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7.2. Disk Emission Measure and Luminosity 

Key signatures of the presence of dense disks include the emission in Balmer hues and 
the IR continuum excess over that from the stellar photosphere. The relation between these 
emission properties, the disk density structure, and the continuum polarization has been 
discussed by various authors (e.g. Waters & Marlborough 1992; Bjorkman & CassineUi 
1990, Yudin 2001, Wood, Bjorkman, & Bjorkman 1997). Here, as another check of our 
derived density distribution, we estimate the line and continuum emission that is expected. 
Both the Balmer lines and the IR continuum emission arise collisionally, so at moderate 
optical depths their luminosities L take the form L = Anj x EM, where EM — n^dV is 
the emission measure of source volume V, and j is the appropriate emission coefficient in 
[ergs cm^ s~^] for each process, which we assume to be constant through the source volume V 
since the cooled post-shock disk is taken to be isothermal. We assume the disk has a roughly 
constant thickness H at all distances from Xinner to Xouter, and use the density distribution 
in Equation (30) for Pb{x), to find 



f-^outer 

EM = 27rR^H I nlxdx (49) 



•dinner 
32 u r „ 2 



2TiR'HlLnj,/ (50) 
II X EM, (51) 



where II is given by 



^L^J^i Pi -dinner ) -^outer) 



/•^^outer / 1 \ 

/ x-2''+Ml--J dx (52) 

'dinner ^ '^Z 



The integral II can be expressed in terms of incomplete beta functions. For integer h and 
half integer (3 it can be evaluated analytically and we find it insensitive to (3 and to XAif ^ 1- 
For 6 = 3 its value is in the range = (1 — 5) x 10~^ for 2 > Xmner > 1- The scale emission 



measure EM^ is 



EMr. 



1 H \ Mv, 



= 2.3 X lO^^cm-2 — 
R 



(53) 
(54) 

(55) 
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If we use an approximate value for the emissivity, 47rjHa = 2.4 x 10 ^^ergs cm^ s ^ (Oster- 
brock 1989) we find the hne luminosity, L(Hq;) = 'i'^jua x EM to be 



L(Ha) = 6.3 X lO'^^erg s 



R 



(56) 



Also of interest is the comparison of the line luminosity L(Ha) with the stellar continuum 
Lx at the Ha line wavelength. For the purposes here the stellar continuum luminosity can be 
approximated by the Rayleigh- Jeans distribution B^^ {T). Thus the monochromatic stellar 
continuum luminosity at Ha is 



T* ~ 



Sw^R^ckT 
A4 



(57) 



This yields the optically thin hne equivalent width Wx — L(lia) / L\{T) using Equations (57) 
and (56) to be 



647r3 c R R^{kTf 



Wx 



(58) 



giving, for II ~ 0.05 



„ H f M 

Wx = i.sA — \ 3-^1 



(59) 



For larger M stars this yields a huge value for Wx since we have ignored attenuation. We 
roughly correct for this by evaluating 



w;^ = Wxe-^ 

and use an approximate disk electron scattering optical depth 



(60) 



T UoRcTt ~ 0.1 



M.gVs 
R12T4 



(61) 



Results for r, VF( versus spectral type Tes are given in Table 2 and shown in Figure 12. 
For stars with higher T^s the high disk optical depth can suppress while for the lower 
Teffthe low density results in small Wx even before optical depth effects are taken into account. 
According to this model, only a narrow range of T^s around that of early B stars should 
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show an easily detectable Ha line. Observationally, Jaschek & Jaschek (1983) found that the 
frequency of Be stars as a function of spectral type peaks at spectral class B2, and a similar 
result is seen in the IRAS survey of B stars by Cote & Waters (1997). Equivalent widths of Ha 
range from a few A to about 30 A (Cote ct al. 1996; Cote & Waters 1997). The Magnetically 
Torqued Disk model predicts roughly the right emission line properties, although clearly a 
better treatment of the radiation transfer is required. Also the Be phenomenon extends to 
late B stars, so a reason for the wide range of spectral classes that contain Be stars needs to 
be explored in the context of our model. We now consider the implications of this model for 
the qualitatively different phenomenon, the variability of the Be stars on a variety of time 
scales. 



7.3. Time Vciriable Phenomena 

Undoubtedly many of the chief clues to the dynamical structure of Be star disks come 
from the complex time variability of these objects. Be stars vary on a wide range of charac- 
teristic times from "activity" on scales of hours (Balona 2000) to quasi-periodic variability 
on years to decades (Doazan 1982). Also of interest is the phase transition time from non- 
emission line to emission line phases that indicate disks can vanish and re-appear on time 
intervals of hundreds of days, (Hubert-Delplace et al. 1982 and Hubert 1981; Dachs 1987). 
The connection between our Magnetically Torqued Disk (MTD) model and these time scales 
is discussed below in § 7.3.1. In addition to changes in overall line strength, another interest- 
ing source of variability in these objects is the cyclic change in the strengths of the violet and 
red peaks of the Ha emission lines, known as V/R variation. This phenomenon comprises a 
shifting asymmetry of emission line peaks from V to R over a timescale of order a decade 
(Telting et al. 1994; Telting 2000). Ideas for understanding the V/R variations within the 
context of the MTD model are discussed in § 7.3.2 



7.3.1. Time Development of Disk Thickness 

Our interest has been in finding the field strength required to have matter transmitted 
from the star to the relatively dense disk inferred from observations. Later we mention 
some ideas for modifying our basic model to produce a gas dominated disk that seems to be 
required for the one armed instability to operate. 

As can be seen in Equation (59), we can attribute time variability of Ha, to the time 
dependent increase in the height of the disk, H{t), so we need an expression relating the 
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time derivative dH/ dt to stellar parameters. We first note that the mass inside a disk area 
element A{d) at any time is p{d)A{d)H . Both p{d) and A{d) should not change appreciably 
in our rigid field model. If we equate the change in this mass to the mass infiux from the 
magnetically channeled wind, we obtain 

pA{d)'^ = Tu{d)A{d) (62) 
where p{d) = Tu{,d)vw{d) / c\ and we get 

^ = ^ = 0.52 X 10^2 cm/year — (63) 

dt Vy, vs 

Let us consider the time for the star to change from the non emission phase to an emission 
phase Be star. For this let us assume the lines are noticeable when the equivalent width is 
about 1 A. Therefore, setting Wx = 1.0 in Equation (59), we obtain from H{t) = {dH/dt)xti, 
the start-up time. 



h = 1.5 years . ^ (64) 

For the values associated with a B2 star, and with T4 = O.STeff this gives a reasonable value 
of about 4 months. For later type B stars, the value increases significantly (100 years in 
the case of a B5V stars and 10^ years at B9V). These long times for the late B stars is 
mostly because the mass loss rates we assumed for these stars are much lower, and hence ti 
in Equation (64) is longer. Although there is a known a tendency for later B stars to change 
on a slower time scale (Hubcrt-Dclplacc ct al. 1981), wc suggest that a better agreement 
of our model with the observations of the late B stars would occur if we accounted for the 
extreme version of FMR theory that was called "Centrifugal Magnetic Rotator" theory (L 
& C, §9.7.1 and Figure 9.7). In this class of magnetic rotator winds the mass loss rate is 
enhanced because the density scale height in the subsonic and solid body portion of the 
wind is increased by centrifugal forces and M increases exponentially with So^ (L & C, Eq 
9.90). The late Be stars appear to have a larger <So on average than do the early Be stars 
(Hubert-Delplace et al. 1981), and would be the ones for which centrifugal magnetic rotator 
effects are likely to be important. If late Be stars have mass loss rates in excess of those 
predicted by pure line driving, our estimate for the disk start-up timescale could be brought 
to an reasonable range. 

In our model, the height of the disk would increase linearly with time, so a growth-time 
scale can be defined as tc = 1/ (dlnH/dt). Of particular interest in this regard is the fill- up 



time tfiii. This is the time it takes the incident mass flux to increase the height H of disk 
matter in the tube to its maximal value which is about d (c.f. Donati et al. 2001 discussion). 
Thus 



tfiii ^ -jh = — — = tpiow [ — ] X (65) 



^ \CsJ \cs 

R12VS 

~ 2.1 years — — — x x (66) 

Made without considering detailed mechanisms, this estimate of the disk filling time is much 
longer than the wind flow time {tpiow = -R/^oo)- The disk filling time is about four years in the 
inner disk (x ~ 2), with the outer regions eventually filling in around a decade — tantalizingly 
close to the observed long term Be star disk variabihty time scales. Some discussion of specific 
field disruption mechanisms (reconnection, diffusion) by matter accumulation is to be found 
in Havnes & Goertz (1984), Babel & Montmerle (1997a), and Donati et al. (2001). 



7.3.2. Connections with the V/R Variations 

The interpretation of the V/R variations proposed by Papaloizou et al. (1992) and 
Okazaki (1991) is that the variation is due to the very slow rotation/precession of an az- 
imuthally localized density enhancement in the disk across the line of sight from the ap- 
proaching to the receding part of the disk. The physical model proposed for such a density 
enhancement and drift is that it is a density wave in the disk, somewhat like galactic spiral 
density waves, but attributable to the gravitational field of the non-spherical central star, 
rather than to self-gravitation (negligible in a stellar disk). This model is attractive in pre- 
dicting a density spiral precession period of the right order for typical Keplerian speeds. The 
MTD model discussed here shares some important features with the models of precessing 
one-armed spirals. Most importantly, both models envision a disk in which material spends 
a long time in roughly Keplerian motion (in fact, the observations of V/R variations are 
often cited as evidence of the existence of a near-Keplerian disk) . 

The V/R variability on irregular time scales of average 7 years (Okazaki 1997) has been 
explained as the orbital motion of a one armed spiral pattern in a disk. A model by Lee, 
Saio, & Osaki (1991) has proposed that the matter is driven up from the star by viscosity. 

A one armed pattern in the disk has been explained by the Global Disk Oscillation model of 
Okazaki (1991, 1996) and Papaloizou et al. (1992). The GDO model would operate in quasi- 
Keplerian disks but not in disks in which the material has a short residence time. Matter 
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must be supplied to the disk either continually or by events such as non-radial pulsations 
that occur on a time scale shorter than the viscous time scale as considered by Owocki & 
Cranmer (2001). 

The V/R phenomenon is one that all models must try to address. In our magnetic chan- 
nehng/torquing model, for strong fields a substantial part of the disk would have corotating 
V(j) values and one might expect the field torquing forces to overwhelm any gravitational den- 
sity wave effect, though in the outer disk where the field will bend (attaining a component) 
there might be some such effect. Any concentration of disk density by a stellar surface field 
that is not uniform in (j) would be expected roughly to corotatc with the star, rather than 
over long periods. This is therefore an intriguing open question for future consideration. 

A very interesting idea, proposed by H. Henrichs (2002, private communication), is 
that the stellar field might be time variable. During periods when it is strong our model 
would apply and a dense disk would build up over the fill time we estimated. However as 
it stands the model provides no obvious way to explain V/R variations, at least via density 
wave precession in a Keplerian disk controlled by the gravity of the non-spherical star. If, 
however, the field weakened substantially then the disk would decouple from the field and 
become near Keplerian through viscous effects and capable of supporting such density waves 
until the field returned or the disk decayed - another observed phenomenon the timcscale 
of which requires explanation. This idea makes the testable prediction that V/R variations 
should be absent during disk build up and present only at times of weak field. 

8. Discussion 

The central point of our model is that magnetic torquing and channeling of wind flow 
from the intermediate latitudes on a B star can, for plausible field strengths, create a dense 
disk a few stellar radii in extent in which the velocity is azimuthal and of order the local 
Keplerian speed. The model was motivated by the fact that line profile data (e.g. Hanuschik 
1996) indicate that disk emission lines of Fe II, extend to Av > VoSini, and reach disk 
Keplerian values. Also the sharp absorption features seen in the center of those lines do 
not show evidence for inflow or outflow. The latter of course would be expected in the case 
of the WCD model. The present paper aimed to show that flelds could conflne and torque 
the flow to create a quasi-Keplerian disk of roughly the spatial extent suggested by these 
observations. Future work to test the real viability of the model will address two important 
questions: 

a) Is the azimuthal velocity distribution v^{d) across the radial range of the dense disk 
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predicted by our magnetic model compatible with observed line profiles? This model 
distribution starts as roughly corotational, i.e. ^,^(0?) oc d, and transitions to a 
decreasing with d. This distribution is clearly not truly Keplerian disk v^{d) oc l/\/d. 
However, the typical is comparable with vk and exists over only a rather narrow 
range of rf, unless the field is very strong, which is unlikely on theoretical grounds 
(Maheswaran & Cassinclli 1988, 1992; c.f. Figure 9). Outside the model disk domain, 
i.e. at X > XAif, we predict quantitatively a somewhat flattened WCD region changing 
with increasing d to a WCZ. At even farther distances from the star, the wind will be 
nearly spherical. All of these regions involve strong outflow velocities but with densities 
much below those in the quasi-Keplerian disk. It is important to estimate whether the 
high Vr but low density distributions there arc compatible with line profile data i.e. 
does there exist an outer expanding flattened disk which affects line proflles very little 
but is detectable by radio observations etc. These and other diagnostic predictions 
of the Magnetically Torqued Disk model such as X-ray production will be treated in 
subsequent papers. 

d) Another issue, that is perhaps related to the V/R variations is the winding up of 
the field - something we have ignored so far. In a model with strictly corotational 
V4>{t) = Vot/R out to some point and angular momentum conservation beyond there, 
the field would corotate in the disk region. However, in practice and for our parametric 
model (Equation 7), v^{r) falls behind corotation and the field will be wound up. The 
resulting component will have opposite sign on opposite sides of the disk plane 
and will grow until field reconnection can occur allowing a quasi-steady state with the 
disk matter slipping across the field lines. In such a situation gravitational control of 
the disk may dominate sufficiently for the disk to behave as quasi-Keplerian and allow 
density waves. Indeed, in the light of this, one might argue that our estimate of i?o,Th 
is too conservative since the crucial feature of our model is that the field should deliver 
wind matter to the disk with sufficient angular momentum (from magnetic torquing) 
and with essentially no Vr but only vq so that shock compression can occur without 
any outfiow (in contrast to the WCD model) . This needs working out in detail but if 
only the wind density needs torquing action our field limits would be much reduced. 
Furthermore, in that case the disk pov'^ would be so high compared to B'^/Sn that the 
matter would completely dominate the field once delivered into the dense disk, which 
would consequently be essentially Keplerian and so support density waves under the 
action of stellar gravity. 

In this paper we have followed a used a kinematic description of the azimuthal veloc- 
ity of a disk zone, along with constraints based on magnetic rotator theory to develop an 
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analytic model which we call the Magnetically Torqued Disk (MTD) model, We derived 
the independent parameters for disk formation problem: the spin rate parameter. So, the 
surface field -Bo.Th and the stellar parameters contained in the quantity, Y^,. We assumed 
a rather simple rigid magnetic field picture, but one that allowed us to derive an estimate 
of the field required to produce a torqued disk for main sequence stars ranging from 03 
V to B9 V. Although an unreasonably large field was found to be needed for the early O 
stars, we concluded that at the early B spectral class at which the Be phenomena are most 
prominent, the model provided about the right observational properties, and for reasonable 
surface field strength. Also, we considered several time scales associated with Be stars, and 
found start-up, growth and fill-up time scales that are in the month to decade long range 
that seem appropriate. For the late B stars the start-up time scale and the Ha equivalent 
widths appear to be small and we suggest that perhaps the mass loss rates for these stars 
need to be enhanced by centrifugal magnetic rotator forces. Modifications that are needed 
to address V/R variations and observational consequences at a variety of wavelengths will 
be considered in later papers. 
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Table 1: Minimum fields for torqued disks around stars along the main sequence 
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Table 2: EqTii\"aleiit AA'idtlis of i/,, from torfiued disks in stars along tlie main sequence 
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Fig. 1. — Illustrates the overall structure assumed in our Magnetically Torqued Disk model. 
The star has a dipole-like inner field with axis aligned with the rotation axis of the star. 
Material leaving the star at latitude / is channeled by the field to the equatorial disk at 
r = d{l). The rotation causes the channeling process to transfer angular momentum from 
the star to the disk. This will lead to a distribution of angular speed f,^ versus the equatorial 
distance d from the star. Note that in the outer regions the field lines are stretched and 
drawn out by the flow from the star. 
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Fig. 2. — The assumed parametric distribution of v^{d) versus equatorial distance. These 
are not flow hnes but rather give the azimuthal velocity ^^((i) of material as it reaches the 
equatorial plane having started at latitude sector I. A spin rate of So — 0.5, is used for 
this figure. In the near solid body region near the star the v^j) for increasing I increases 
with distance, d{l) from the star. Note that a depends on the magnetic field, the spin So 
and outflow density in the wind. Curves for flvc stars of different a values are shown (from 
Equation 7). Also shown are the Keplerian vxid) distribution and the parallel but larger 
escape speed distribution VE{d). 
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Fig. 3. — The curves of S{x) = v^{x)/vk{x) for a range of a and for So = 0.5. To contribute 
to a quasi-Keplerian disk the S distribution must reach values of unity and above. There 
is a minimal value of a for this to occur. For higher values of a, S reaches a maximum 
farther out and that maximal point is taken to be the outer boundary of the magnetically 
dominated disk. The line types are identified in Figure 2. 
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Fig. 4. — The minimal and maximal radii dinner and Xouter as functions of a. Equation 10 is 
used to calculate dinner with the indicated values of So, while Xouter = V^a. 
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Fig. 5. — The solid line shows the dependence on the spin parameter So of the minimal value 
of a (= ajh) required for disk formation, as given by Eq. (15). The dashed line illustrates 
the good fit given by the power law expression described by Eq. (16). 
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Fig. 6. — The function r{a) (Eq. (36)) is shown for the values of b and /3 as described by 
the different hne styles in the legend. The case 6 = 3 corresponds to a dipole magnetic field. 
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So 

Fig. 7. — The threshold value of F versus So for 6=3 and for /3 = 1.0, from empirical 
Equation (42) (dashed line) and also from the exact solution of Equations (15) and (35) 
(solid hue). 
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0.2 0.4 0.6 0.8 1.0 

So 

Fig. 8. — The minimal or threshold surface magnetic field for the formation of a disk versus 
»So for which we have used the parameters of a B2V star from Table 1. The different curves 
illustrate the effects of varying the model parameters, b, which gives the radial dependence 
of the magnetic field, and /3, which is the outfiow velocity law parameter. For a case in which 
the field drops more sharply with distance from the star the required surface field needs to 
be larger. Three values of f3 are shown for the 6 = 3 case. The slower the rise in velocity, 
(e.g. (3 — 2) the smaller the field needed to produce the torqued disk, for a given So- 
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Fig. 9. — Shows the hmits derived for the maximal and minimal field strengths for Be stars 
versus So derived from photospheric constraints of Maheswaran & Cassinelli (1992). The 
solid line shows the minimal Bo as a function of So such that a disk can be magnetically 
torqued; for these calculations we use the stellar parameters from Table 1,6 = 3, and /? = 1. 
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Fig. 10. — Two panels that show the radial extent (dinner to Xoutcr) of magnetically torqued 
disks versus the surface magnetic field, B^. The B field is normalized in units of F*, so the 
plot is apphcable to any star, a) In the left panel results for the boundaries are shown for 
three values of the spin rate parameter S^- The upper segments of the curves correspond to 
Xouter and note these merge with the lower segments, the dinner curves, at the threshold point 
for each of the three models. The curves appear very similar in shape although shifted, b) 
To illustrate the homologous nature of the curves, in the right panel is shown the quantities 
on both axes in the left panel are normalized by the threshold values. Notice that the extent 
of the disk (xouter- ^^inner) iucrcascs sharply for magnetic fields just above the threshold value. 
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Fig. 11. — Shows the threshold field for torqued disks for stars along the main sequence. 
The stellar parameters are from Table 1 and b — 3, P — 1 are assumed. The three curves 
correspond to the values of the spin rate So — 0.5, 0.7 and 0.9, two of which are labeled. 
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Fig. 12. — The soHd Hne shows an estimate of the upper limit to the equivalent width of the 
Ha line in A versus T^s, produced in a disk of infinite extent. The dotted line shows the 
reduced value of the equivalent width when optical depths in the disk are accounted for in 
an approximate way. The unmodified Wx indicates unacceptably large values for the early 
type stars, however, for these stars the magnetic field required would also be excessively 
large (i.e. greater than a kilogauss). 



